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1. Introduction 

The behaviour of QCD at non-zero temperature and density is an active area of 
research both from a theoretical as well as an experimental perspective. While per- 
turbation theory can be used at very large temperatures or densities, this is of limited 
use in intermediate regimes such as those that are relevant for heavy ion collisions 
or the description of the interior of neutron stars, where QCD is strongly coupled. 
Lattice methods in turn are hard to apply to highly time-dependent processes, and 
also pose difficulties when dealing with systems with large baryon density because of 
the sign problem. Alternative techniques are therefore called for. The gauge/gravity 
duality provides a non-conventional approach to the study of strongly coupled gauge 
theories, and while the dual of QCD proper is not known (not even in the large- iV^ 
limit), some of the current proposals for holographic duals of gauge theories do at 
least seem to capture several qualitative properties of strongly coupled QCD. 
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One of the key questions related to QCD at non-zero quark chemical potential is 
that of the character of the ground state. Various models suggest that, for particular 
values of the chemical potential and temperature, an instability sets in which results 
in the decay of the homogeneous vacuum to one which breaks translational invariance 
(see for instance the ladder analysis of QCD in the large-A^c limit by Deryagin, 
Grigoriev and Rubakov [1]). In the dual holographic language, chemical potentials 
are typically realised by introducing an electric gauge field configuration on the world- 
volume of a flavour brane, with the non-normalisable component being related to 
the value of the chemical potential. Such electric configurations in the presence of a 
curved background can become unstable, as was emphasised in work of Nakamura, 
Park and Ooguri [2]. These authors have pointed out that a black hole with a 
constant electric field in five-dimensional Maxwell- Chern-Simons theory is unstable 
against decay towards a configuration with both electric and magnetic fields turned 
on. Importantly, this instability occurs only for a non-vanishing value of the spatial 
momentum k, signalling that the new ground state is spatially modulated. The 
lowest energy state has a uniquely specified momentum k^s, which is determined by 
the value of the electric field. The origin of the instability can be traced to the 
presence of a Chern-Simons coupling. 

In the context of holographic models, such behaviour was in fact first found in 
the work of Harvey and Domokos [3], who exhibited an instability in a particular 
bottom-up model with a Chern-Simons coupling and an axial chemical potential. 
More recently, Ooguri and Park [4, 5] have analysed the Sakai-Sugimoto model [6, 7] 
in the deconfined phase, where chiral symmetry is restored. This model contains 
Chern-Simons couplings in the action of fiavour D8-branes. The presence of a non- 
vanishing electric field corresponding to a quark chemical potential then leads to 
an instability towards a phase where baryonic charge density is generated. This 
non-homogeneous ground state contains corresponding baryonic and axial currents 
carrying non-zero momenta kgg. 

This still leaves open the question as to whether a similar instability can occur 
at vanishing temperature. In the aforementioned Sakai-Sugimoto model, it is known 
that even at zero temperature, a sufficiently large isospin chemical potential can lead 
to an instability towards a new ground state in which vector mesons condense [8]. 
This ground state, even though it breaks rotational invariance, is spatially homoge- 
neous. The question remains whether a non-homogeneous ground state is possible in 
this model. In the present paper, we therefore set out to use the holographic model 
of Sakai and Sugimoto to study the properties of a QCD-like theory in the large- A^c 
limit, at non- vanishing values of the axial chemical potential.^ More precisely, we 

^An axial chemical potential is induced by a time-dependent theta angle, and as such is one of 
the motivations to study the Chiral Magnetic Effect [9], where theta angle fluctuations together 
with a magnetic field induce an electric current. The axial potential also plays an important role 
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analyse the ground state in the confining phase where chiral symmetry is broken. 
Ahhough the axial currents are anomalous in real QCD, this anomaly is suppressed 
in the large- iVc limit, and the axial U{1) current is conserved. We will see that this 
potential can trigger an instability to a new non-homogeneous ground state. This is 
to be contrasted with the effect of a baryon chemical potential. The latter would lead 
to condensation of baryons, which are somewhat problematic in the large- A^c limit as 
they become infinitely heavy, and are in fact singular objects in the Sakai-Sugimoto 
model. ^ Our solutions are regular and have finite energy.^ 

Our results show that even the low-temperature phase of the Sakai-Sugimoto 
is unstable against decay into a non-homogeneous state. These results can be con- 
trasted with recent work [12] on the efi^ects of the axial chemical potential in the 
PLSMg model, where no evidence for an axial potential induced phase transition at 
zero temperature was found. The spatial modulation of the ground state we find 
is characterised by the momentum, which is interestingly, to good approximation 
independent of the axial density (or chemical potential) . We also find that the axial 
density of the new ground state is substantially larger for a given chemical potential 
than in the homogeneous phase. The instabilities towards new ground state occur 
both in the Maxwell-Chern-Simons truncation of the flavour D8-brane action as well 
as in the full DBI system, provided the Chern-Simons coupling of the latter is suf- 
ficiently large. We do not find any evidence for condensation into other states as 
the chemical potential is increased further. While new non-homogeneous configura- 
tions appear, they do not correspond to new global energy minima, but only to new 
meta-stable vacua, as they have larger free energy than the true ground state. 

2. Chemical potentials in the Sakai-Sugimoto model 

2.1 Review of the model 

In this section we briefly review the basics of Sakai-Sugimoto model [6, 7] in order to 
set the scene and to introduce the conventions we will be using in this paper. The 
Sakai-Sugimoto model at low temperature is derived by considering the decoupling 
limit of a large number Nc of D4-branes compactified on a circle of radius R spanned 
by the coordinate r. Anti-periodic boundary conditions are imposed on the fermions. 
The system furthermore includes Nf fiavour D8-branes positioned at r = and 
Nf anti-D8-branes positioned at r = L. The dual gauge theory is a maximally 

in the holographic realisation of chiral magnetic spirals [10]. In the present paper, however, we do 
not turn on any additional magnetic field but simply study the consequence of the axial chemical 
potential by itself. 

^Thc authors of [11] found a dynamical instability treating the baryonic chemical potential as a 
source of cusps in the flavoiu D8-branes. 

^Because the solutions originate from the Chern-Simons term, they are of course still such that 
the fields scale as ~ A, i.e. their scale is set by the string length. 
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supersymmetric SU{Nc) gauge theory in 4+1 dimensions, which is compactified on a 
circle with anti-periodic boundary conditions for the adjoint fermions, coupled to Nf 
left-handed fermions in the fundamental representation of SU (Nc) localised at r = 
and Nf right-handed fermions in the fundamental representation localised at r = L. 

In more detail, the background generated by the stack of D4-branes has a metric, 
Ramond-Ramond four-form and string coupling given by 

TT3/2 7-J-3/2 p3/2 ]tt2 

(2.1) 

where 64 is the volume form on 5*^ and V^i its volume. The index /i = 0, 1, 2, 3 defines 
the directions of the dual gauge theory, and r is a coordinate on the circle so that 
r = r + 5t. The curvature radius is given by i? = {TcggNcY^^ls where is the string 
coupling and Ig = \fa' is the string length. The period of the compact coordinate 
introduces a four-dimensional mass scale 

and the four dimensional 't Hooft constant is related to the string coupling by 

A = gl^N^ = 27rMj,^gsN, h . (2.3) 

For our purposes it will also be convenient to introduce use pair of dimensionless 
coordinates y and z defined by the relations 

U = U,^{l + y' + z'y^' = U,,Kl/J , r = ^arctan (^^^ . (2.4) 

In terms of these coordinates the metric takes the form 

+ {y' + ^'Kl/!W + 2y^(l - Kll!)dydz] + R'l'U'J^KlUdnl . (2.5) 

In this background one introduces a probe D8-brane and D8-brane which fill 
out the 3 + 1 dimensional space-time and the four-sphere, and trace out a curve in 
the remaining two directions {U.r). The latter coordinates span a two dimensional 
cigar-like subspace. The shape of the curve is obtained by solving the equation of 
motion of the probe branes, and is given by a one parameter family, which is specified 
by the asymptotic separation between the D8 and D8-brane. In this paper we will be 
interested in the so-called antipodal brane embedding corresponding to the maximal. 
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antipodal separation of the branes in the circular direction r. In this case, the branes 
extend all the way to the tip of the cigar where they join in a smooth way. For this 
special case the metric on the brane world-volume takes the form 

9 

= GA/ivdx^^dx^ , 

(2.6) 

where = 1 + . 

The action for the U (Nf) gauge fields living on the D8-D8-branes is given by the 
DBI action which in the non-abelian limit is not fully known, but at the quadratic 
order reduces just to the the Yang-Mills action in a curved background, corresponding 
to the curved world-volume of the brane, 

^YM = -/U8(W)' j dVd^d^fie-'^v^^Tr (J^^'^'J^mn) (2.7) 

where J'mn = Qm-A-n — O^Am + i[AM,AN], and the index m = {n,z,a). The 
coordinate z spans the holographic direction and a labels the directions on S"^; indices 
are raised or lowered using the metric (2.6). In order to avoid non-singlet 5*0(5) 
states associated to the S^, in what follows we will set all the excitations along the 
sphere directions to zero, Aa = and assume that A^ and Az do not depend on the 
S*^ coordinates. Then we can integrate out the S*^ coordinates and the Yang-Mills 
action reduces to 



5'ym = — K / d xdzTr 

where 



(2.8) 



K = = . (^.y) 

1087r2 2167r3 ^ ^ 

Note that the Yang-Mills action can now be written in terms of an effective five- 
dimensional metric 

dsjs) = gmndx^dx^ = Ml^Kl'^ T^^Ax^dx^ + ^,"2/3^^2 ^ (3. 10) 

in terms of which it takes the form 



SYM = -'^j d^xdz^ Tr [^™" J-^J , (2.11) 

where the indices are now raised or lowered using the metric Qmn- 

The effective action for the D8 brane includes also a Chern-Simons term which 
is given by 

Scs = f^8^^;r^ A = a , a = a—^ . (2.12) 



3! .Ids ^ ./ ^ ' 24n 
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Here d = 1 in string theory but we will keep it as a separate parameter for later 
use^. Again, is the Ramond-Ramond four-form, proportional to the volume form 
of S'^ (see (2.1)), and the form is the usual five-dimensional Chern-Simons form, 

CU5 = Tr (^A:F^ - '-A^7 - ^^A^ . (2.13) 

It is useful for our purposes to decompose the U{Nf) gauge field into the SU{Nf) 
part A and the U{1) part A as 

A = A+ ^^A (2.14) 

For the case Nf = 2, we can decompose the U{2) symmetry into the baryon U{1) 
and isospin SU{2), and the Chern-Simons term decomposes as 

Scs = (y j {^ATiJ^^ + ^^/2 + (total derivatives)^ , (2.15) 

We see that in the case of two D-branes, the Chern-Simons coupling has two terms, 
one involving only U{1) fields and the other one which couples the U{\) and the 
SU{2) fields. In the present paper we will focus on the effect of the U{1) part. 

2.2 Holographic realisation of chemical potentials 

Holographic models encode global symmetries of the dual gauge theory in the form 
of gauge symmetries in the bulk theory, and these relations hold both for the closed 
as well as the open sectors of the string theory. So in particular, in the D8-D8 
system, there are two independent gauge theories living near the two boundaries of 
the flavour brane system, giving rise to U{Nf)L and U{Nf)fi gauge symmetries, and 
these correspond to the global U{Nf) x U{Nf)ji of the dual gauge theory. 

In the low-temperature phase of the Sakai-Sugimoto model that we are interested 
in, the two branes are connected in the interior of the bulk space, and thus the gauge 
fields Am and A^ are limits of a single gauge field living on the two connected branes. 
Therefore, one cannot independently perform gauge transformations on these two 
gauge fields, but is constrained to gauge transformations which, as z ^ ±oo, act in 
a related way. Specifically, since near the boundary a large bulk gauge transformation 
acts as Al/r — )■ gL/RAL/Rg^^jj^, then clearly if gL = gR. any state (and in particular 
the trivial vacuum ^ = 0) is invariant under the vectorial transformations, i.e. those 
transformations for which gi = gR- Hence it means that the vector-like symmetry is 
unbroken in this model. 

The correspondence between the bulk gauge field and the source and global 
symmetry current of the dual gauge theory is encoded in the asymptotic behaviour 



*For comparison, our a corresponds to 4/3 times a used in [5]. 
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of the former. More precisely, the bulk gauge field a;^) behaves, near the 

boundary and in the = gauge, as 

z) ^ a„(x'^) (l + 0{z~^/^)) + pn{xn\[l + 0{z-^'^)) . (2.16) 

Here a„(x^) describes non-normalisable behaviour of the field, and is interpreted as 
a source in a dual field theory action of the form 

(2.7) 

while Pn{x^) is proportional to the expectation value of the current Jn{x^). Therefore, 
adding a chemical potential to the field theory corresponds to adding a source for J°, 
which implies the boundary condition for the holographic gauge field An{x) = /i5„o- 

For the Sakai-Sugimoto model, the bulk field Am living on the D8-branes has two 
asymptotic regions, corresponding to each brane, and hence there are two indepen- 
dent chemical potentials and pr which can be separately turned on. Instead of 
left and right chemical potentials one often introduces vectorial and axial potentials, 
defined respectively as pv = + Pr) and fiA = ^{pl — Pr)- The vectorial and 
axial chemical potentials for the U{1) subgroup of the U{2) gauge group on the two 
D8-branes correspond to the baryonic and axial chemical potential in the dual gauge 
theory, while the non-abelian SU{2) chemical potentials are mapped to the vectorial 
and axial isospin potentials. 

In this paper we will be interested only in axial chemical potentials. Our compu- 
tations will mainly be concerned with the U{1) subgroup, since the instability which 
leads towards the spatially modulated phase in the U{1) subgroup is present in the 
same form in the U{1) subgroup of an isospin SU{2) group, leading to the instability 
of homogeneous configuration for the isospin chemical potential as well. 

Our starting point is the homogeneous solution [6] in the Az = gauge. For the 
configuration to describe an axial chemical potential one needs to turn on only the 
Ao{z) component, and look for odd, non-normalisable solutions to the equations of 
motion (3.2). It is not hard to see that with this simple ansatz the Chern-Simons 
contribution to the equations of motion vanishes, and the solution to the equation 
of motion is given by 

2 

Ao{z) = —fiASi.Tcta.nz . (2.18) 

71 

Just as in the chiral Lagrangian, there is no potential generated on moduli space [8] , 
and the result of the axial potential is simply a non-vanishing axial density. The 
value of this density is easily read off using (2.16) to be pA = ^pa/t^- This same 
solution can also be trivially embedded in the U{1) subgroup of the SU(2) isospin 
group. In this case however, as explained in [8], the odd Ao{z) configuration in the 
Az = gauge is equivalent (by a global SU (2) rotation) to a system with a vectorial 
chemical potential and a non-trivial pion condensate. 
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3. The spatially modulated phase 

We have seen in the previous section that turning on the axial chemical potential 
in Sakai-Sugimoto model requires finding a non-normalisable spatially homogeneous 
solution for the U{1) field Aq, which is an odd function in the holographic direction. 
The asymptotic values of correspond to the values of the chemical potentials in 
the dual gauge theory. 

As mentioned in the introduction, the work of Domokos and Harvey [3] and 
Nakamura, Ooguri and Park [2, 4] has shown that in Maxwell theory with a Chern- 
Simons coupling turned on, a constant electric field is unstable against decay to a 
spatially modulated phase, in which a magnetic field transverse to the direction of 
the initial electric field is switched on. Since the configuration with non-vanishing 
axial chemical potential (2.18) amounts to having a non- vanishing electric field in 
the bulk, we expect that an instability similar to that found by the aforementioned 
authors should be present here. In the present section we will analyse the Sakai- 
Sugimoto model at zero temperature with non-zero axial chemical potential, show 
that for sufficiently large values of that potential a new ground state appears, and 
determine that it has lower energy than the homogeneous vacuum. 

3.1 Equations of motion 

Let us begin by writing the effective 5-dimensional action for the case of a single 
fiavour Nf = 1, obtained after we integrate out the fields on the four sphere S^, for 
the case of the single D8-D8 brane, 

^YM + 5cs = -| y d^xdzv^^"^"^-^^ + I e^™ j d^xdzAeJ'mn^Fp, . (3.1) 

The equations of motion following from this action read 

^ Kj K 

(3.2) 

Here we have split the five- dimensional bulk indices as m = {i = (1, 2, 3), 0, 2;) and 
we have assumed that the general diagonal five-dimensional effective metric gmn 
depends only on the z coordinate and has symmetry gn = g22 = gss = gxx (valid for 
the Sakai-Sugimoto system). 

It is important to note that the variation of the action (3.1) leads to a boundary 
term, which, after one uses the Bianchi identity dmJ^pq + dqj^mp + dpj^qm = 0, reads 

55bdy = - / d^xdz dm [SAe {2k^ J^"^^ + ae^™A-7>c)] • (3.3) 
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Therefore, when 6Ai does not vanish, one has to add a boundary term to the action 
in order to ensure that there are no boundary contributions to the equations of 
motion (see e.g. [13] for a discussion in related context). This is important when 
we consider the system at constant density. In the appendix we discuss this issue in 
more detail, and give the corresponding boundary term for the antipodal embedding 
of D8-D8 branes. 

3.2 The abelian ansatz 

To describe a spatially modulated phase, our starting point is the ansatz 

Ao = fiz,x), A = Aiz,x), (3.4) 

— * 

where A are in the spatial directions of the boundary. We will work in the gauge Az = 
0. Furthermore, we only want to turn on the chemical potential for the axial currents 
and insist that the spatially modulated phase is spontaneously generated (i.e. without 
introducing any kind of sources). Hence we impose normalisable boundary conditions 
on the fields A, while similar to the homogeneous case we impose that is non- 
normalisable, and tends to opposite constants at the two ends of the D8-brane system. 



A{z = ±oo) = 0, Ao{z = ±oo) = ±-iJ, A- (3.5) 

TT 

For this ansatz (3.4) the full non-linear equations of motion (3.2) become^ 



a 



-9 9^'g''dz{V-A)-3-£-B 







a 



dz 



-gg'^g'^'^dzA 



^g^^g--{W-£)-3-idzA)-B = 0, 

K 

^g{gn^V xB + 3- UdJ)B- dj x £ 



(3.6) 
(3.7) 
0. (3.8) 



In these equations, S = J^oiX^ = —V/ and B = \eijkJ^jkX^ = V x ^ are the (bulk) 
transverse electric and magnetic field associated to / and A. 
A simple ansatz that solves these equations is to impose 



£ = 0, B = kA. 



(3.9) 



where k = ±|A;| and k is the spatial momentum. This way the equation (3.6) is 
automatically satisfied and the equation (3.7) can be integrated. The gauge field 
consistent with this ansatz can be written as 



^0 = f{z) , A = f]i{z) cos{k ■ x) + f]2{z) sm{k ■ x) 



(3.10) 



^In order to compare our normalisation of the Chcrn-Simons coupling a with [5], note that 
equation (6) of that paper contains a typo: the coefficient a should read a/6. 
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where fji, fj2 are polarisation vectors satisfying 



fji X k = kf]2 , k X f]2 = kfji 



(3.11) 



For simplicity we choose k = k xi so that the polarisation vectors reduce to ffi{z) = 
h{z)x2, '>]2{z) = —h{z)xs where h{z) is a function that vanishes at z ^ ±00. The 
equations (3.7) and (3.8) take the form 



3 a 
2 k 



Ml^KA [KAh] - Ml^Kf'Ph + 3-kh 



K 



kh' 



P 



3 a 
2 k 



kh' 



0, 



(3.12) 



(3.13) 



where p is an integration constant, = 1 + z"^ and M^k is defined in (2.2). Equa- 
tions (3.12) and (3.13) form the system that determines the new vacuum. In order 
to solve numerically this system it is convenient to introduce dimensionless variables 
/, h, k and p defined by 



f = \M^^f h = XM^^h k = M^J , p = \Ml^p, 
with A = A/(277r). Then equations (3.12) and (3.13) can be written as 



KAf = P 



(3.14) 



(3.15) 



KA 



K,dM 



a 



kp - i^2/3p 



h-—k''h' 
2 



0, 



(3.16) 



To determine the physical meaning of the integration constant p, we recall that 
the charge density of the dual axial current is computed using the AdS/CFT 
dictionary as 



PA = {Jo)± ~ lim z^d.Ao + lim z^d.Ao ■ (3.17) 

2— >oo z— >— 00 



From the asymptotic behaviour h{z) ~ ao/z + ■ ■ ■ , which can be obtained by per- 
turbatively solving (3.16), one deduces that the function h{z) does not contribute 
to the charge so p is indeed the charge density of the dual axial current (up to a 
constant). However, this does not mean that the charge density is the same as in 
the homogeneous case. The reason is that the chemical potential is determined via 



lim ^0 ~ 1™ ^0) 



dzd;,Ao(z) . 



(3.18) 



Therefore, we see that the constant ao does contribute to the value of the chemical 
potential, and hence modifies it in a nontrivial way so as to give a new relation 
between the chemical potential pa and axial charge pa- We will see this explicitly 
in the next section (see figure 3). 
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Figure 1: Result of the stability analysis in the Maxwell-Chern-Simons system. The 
curve shows the values of p and k for which an instability sets in. From this analysis one 
finds pcrit ~ 2.35, and the corresponding mode is vectorial. The dashed and dotted curves 
correspond to additional instabilities of axial and vectorial nature respectively, which set 
in at higher values of the density. 

3.3 Analysis of the ground state 

In order to find the spatiaUy modulated phase (3.10) we need to solve equation (3.13) 
for h{z) and then use this solution when determining f{z) using equation (3.12). 
Equation (3.13) can be solved numerically using a shooting method. 

Before presenting our numerical solutions, let us note that equations (3.13) 
and (3.12) have three free parameters: the Chern-Simons coupling a defined in (2.12), 
the wave number k which sets the frequency of the spatial modulation, and the 
charge density of the dual current p which appears as an integration constant in 
equation (3.12). The wave number k will be kept free, i.e. we solve the equation for 
various values of this parameter. In string theory the constant a has a fixed value, 
but one might consider a more phenomenological attitude in which a is allowed to 
take on any value. In this section we will set a to its string theory value, so that 
a = 1. The reason is that in Maxwell-Chern-Simons theory, which we discuss in 
this section, the actual value of coupling is not really relevant: for any a, a non- 
homogeneous state can be made to appear since one can always make the chemical 
potential large enough. The situation is however different when dealing with the DBI 
action, which we discuss in §4. 

Because we do not want to introduce any external sources in the boundary theory, 
we are looking for normalisable solutions of the equation (3.13). This means that. 
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Figure 2: The charge density of the non-homogeneous solution, as weh as the energy 
difference with respect to the homogeneous vacuum, as a function of the momentum. These 
plots show the range of momentum values {k^m, ^max) for which non-homogeneous solutions 
exist. Both plots are made for a charge density p = 3.9. 

even though (3.13) is a second order equation, there is only one undetermined param- 
eter, ao, which governs the behaviour of the function h{z) near infinity. By expanding 
the equation (3.13) around ±oo it is easy to see that h behaves as h{z) ~ ao/z + • • • . 
Using the AdS/CFT dictionary we see that the constant ao is related to the expec- 
tation value of the dual current in the directions ffi and 172 on the boundary, 



If the current is vectorial, then the same values are obtained by taking the limit 
z — ±00, otherwise it is axial. 

Next, we use a shooting method to numerically solve (3.13) for various values of 
charge densities p and momenta k. The first observation is that there is a critical 
value of the charge density (and corresponding potential) so that if the density is 
larger than pcrit, non- homogeneous solutions will exist. This critical value is best 
obtained using a linear analysis, looking for unstable modes in the homogeneous 
vacuum. A fluctuation with frequency u satisfies the linear truncation of (3.13) 
with an additional term proportional to u added. A marginally unstable mode is 
then obtained by solving for modes which have a; = 0, i.e. modes satisfying the 
linear truncation of (3.13). Such modes occur for one-dimensional subspaces of the 
p — k plane. Figure 1 depicts the result of this analysis, which shows that there are 
various branches of unstable modes. More specifically, we have found that for any 
P > Pcrit ~ 2.35, instabilities exist. 




(3.19) 
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Figure 3: The chemical potential as a function of the charge density for the non- 
homogeneous ground state. The thick dashed line describes the higher-energy, axial branch 
of the solutions, which are meta-stable. The thin diagonal dashed line shows the relation 
fx A = 7r/5^/2 for the homogeneous vacuum. 



Corresponding to these instabilities are new ground states. For any momentum 
in the range [fcminlp), ^max(p)], where k^i^ and fcmax are the two values of A; on a 
line of instability at fixed p, we find non-homogeneous solutions of the non-linear 
equations. See figure (2) for an example at p = 3.9. However, among all solutions for 
a given p, there is only one momentum k = kgs{p) for which the solution has minimal 
energy. Hence this value specifies the unique ground state for a given charge density 
or chemical potential. We should emphasise that when we determine the ground 
state, we have compared the homogeneous and non-homogeneous solutions at fixed 
charge density p. In other words we are working in the canonical ensemble. One 
could also work in the grand canonical ensemble, comparing solutions at the same 
value of the chemical potential p^. The results which are obtained in that case are 
qualitatively the same as in the canonical ensemble, so we choose to present only the 
canonical ensemble results here. 

The solutions describing ground state are uniquely specified by the value of the 
charge density p. In other words, the parameters (ao; kgg, Pa) specifying the ground 
state are all functions of the charge density p. The dependence of the chemical 
potential on the charge density is plotted in figure 3, which also shows the relation 
between these two parameters in the homogeneous solution. We see that the relation 
is again (almost) linear in the non-homogeneous case, but for a given value of the 
charge density the chemical potential in the homogeneous case is larger than in the 
non-homogeneous case. The dependence of the momentum A;gs on the density p is 
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Figure 4: The momentum of the spatial modulation as a function of the charge density 
for the non-homogeneous ground state. The dashed line describes the higher energy, axial 
branch. Observe that the momentum is practically independent of the density. 



shown in figure 4. Interestingly, we see that the momentum is almost independent 
of the charge density. It would be interesting to understand better the reason for 
this behaviour. To see if this is a consequence of the Maxwell approximation of the 
DBl action, we will examine this relation again in the nonlinear DBI case in the 
next section. Finally, the dependence of the constant on the density p is shown 
in figure 5. As is visible from equation (3.19) this constant is equal to the amplitude 
of the current components J^^ and Jr^^. 

The difference between the free energies of the homogeneous and non-homo- 
geneous configurations (see equation (6.5)) for fixed value of the charge density p 
is displayed in figure 6. As required, we see that when p > pcrit, the free energy 
of the non-homogeneous configuration is lower than the one of the homogeneous 
configuration, indicating that the system will settle in the new ground state. We 
also see that for large enough charge density, p > 5.4, a new branch of solutions 
appears. While the first non-homogeneous branch is given by an even function h{z), 
thus describing a vectorial current, the next branch is an odd function, describing an 
axial current. As the chemical potential is increased further, yet more branches occur, 
alternating between vectorial and axial. However, as is manifest from figure 6, all 
higher branches have higher values of the free energy, and thus correspond to excited 
states, not ground states. 

Finally we should note that while the ground state is characterised by vectorial 
spatial components rji and 772 of the current J, the charge density is always axial. 
This is a consequence of the fact that the integrand in (3.12) is an even function for 
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Figure 5: The amplitude of the expectation value of the currents {Jrn) or (Jjja) ^ 
function of the density. The dashed curve is again the unstable axial branch. 
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Figure 6: The difference between the free energies of the homogeneous and non- 
homogeneous configurations as a function of the charge density. The dashed line depicts 
the axial ground state, which clearly has higher energy and is therefore unstable. 



any function h{z). 

In conclusion, we see that the Sakai-Sugimoto model predicts that, even at zero 
temperature, a non-trivial condensate of vector mesons forms for sufficiently large 
value of the axial chemical potential fiA- This condensate is non-homogeneous, and its 
spatial modulation has a momentum vector which is almost constant as a function of 
the chemical potential. These results complement those in deconfined phase obtained 
by [5]. 
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4. Spatially modulated phase for the DBI action 



In the Maxwell-Chern-Simons theory, an instabihty towards a new ground state of 
the type found in the previous section will always occur, independent of the strength 
of the Chern-Simons coupling. This is simply because the electric field (and hence 
the axial chemical potential (ia) can be made arbitrarily large. However, given the 
fact that the Chern-Simons term is of order with respect to the Maxwell term, 
the fields in this new ground state are necessarily of order A. The solution found in 
the previous section can therefore not be seen as an approximate solution to the full 
DBI system, since the truncation that leads to the Maxwell action assumes that A is 
small. Re-analysing the computation for the DBI action may thus yield qualitatively 
new results.^ 

Furthermore, in the full DBI system, there is an upper bound on the electric 
field, which may prevent the instability from setting in [5]. A final reason to repeat 
the previous analysis for the DBI-Chern-Simons action is to examine the question 
whether the non-linearities of the DBI action modify the values of the ground state 
momentum kgs, perhaps making it dependent on the charge density. 

Our starting point is the DBI action for the D8-D8 system, after integration over 
the S'^ (on which none of the fields which are turned on depend). 



5i 



DBI 



where 
with 



E = det{Emn) , Emn = 9mn + P{z)J^mn , (4.2) 



7(.) = Vs.^^,e-^a{zf'%{zf = ^K;V3 , 



^ ^ a(z) 2A " 



(4.3) 



Here Qmn is the effective 5-d metric defined in (2.10) and a{z) and b{z) are given by 
a{z) = ^M,,R^K;'/' , b{z) = '^M.^R'K'J' . (4.4) 
This action is coupled to the Chern-Simons term 

^ ^imnpq I i4 



^CS = ^ e''""^^ J d^xdz AiJ^.nnJ'pg • (4.5) 

The coupling a is defined to be a = aNc/ {247r'^) with a fixed in string theory to 
a = 1. However, since we are interested in investigating the nature of the instability 
we will relax this condition and consider an arbitrary a. 

^Note that because the momentum scale k of the spatial modulation is set by Mkk (as follows 
from (3.14)), higher-derivative corrections can be ignored when IsMkk is sufficiently small. 
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The action variation leads to the DBI-CS equations 



(4.6) 



and the boundary term 

^'S'bdy = — I d'^xdz dr, 



(4.7) 



Here we have introduced the antisymmetric tensor E^^'"^^ = E^"^ — E^^ where E^^ 
is defined by E^'^Emn = 5^ . 

For the ansatz (3.10), the nonzero components of can be written as 



2(3g^^g^-d,A, 

l + g^-g^op^idjy + g^^g--id,hy' 
2P{g^^ydiA, 



1 + (c/^^)2/32fc2/i2 ■ 



(4.8) 



Using these results the DBI-CS equations can be written in a form which puts all 
modifications due to the DEI action in a single function Q{z): the equation (3.12) 
becomes 

Q{z)K,dJ = p-^kh\ (4.9) 
while the analogue of equation (3.13) for h{z) now reads 



Q{z)KA [Q{z)KAh 
The function Q{z) itself is given by 



+ 



a 



kp - Kl'^k^ 



h-^k^h' = 0, 
2 



(4.10) 




1 + lK'J\dM' 



(4.11) 



and we used again the dimensionless variables introduced in (3.14). 



4.1 Stability analysis 

Before solving the non-linear equation (4.10), it is useful to analyse again possible 
instabilities in the fiuctuation spectrum, in particular in order to see how they depend 
on the Chern-Simons coupling a. Similar to the discussion for the Maxwell case, 
we can find unstable modes by solving for h{z) in the linearised version of (4.10). 
Results for two values of a are given in figure 7. One observes that the number of 
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Figure 7: The relation between p and k for the unstable mode of the DBI system, for 
two values of the Chern-Simons coupling: a = 1 for the left panel and a = 2 for the right 
panel. 
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Figure 8: The relation between the critical value of the density and the Chern-Simons 
parameter a. The critical density diverges as a — t- 1/3. 

solutions which survive from the Maxwell truncation is determined by the Chern- 
Simons coupling, and larger values typically make the curves shift to the bottom left 
of the graph. For a = 1 (the value used in the Maxwell case), the critical density is 
now found to be pcrit ~ 3.46, i.e. larger than in the Maxwell case. 

From this instability analysis one can also find the minimum value of a for which 
an instability is possible at all. Figure 8 shows the critical density pcrit as a function 
of a. This plot agrees with a divergence as d — )■ 1/3, just as in [5], indicating that 
there is no instability in the system when a takes on smaller values. For such values, 
the non-linearities of the DBI action eliminate the instabilities visible in the Maxwell 
truncation. 
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Figure 9: Summary of the DBI analysis, with red (dark) curves depicting the DBI results 
and the Maxwell result is given in orange (light) for comparison. The left panel shows 
the expectation value of the current (J) as a function of the density. The critical density 
has now shifted to pcrit ~ 3.46. The right panel shows that the momentum vector is still 
practically independent of p, but has also changed value. 

4.2 Ground state: Maxwell versus DBI 

Finally, we present the analysis of the non-homogeneous ground state in the DBI 
case. The logic of obtaining it is again the same as in the Maxwell case, so we will 
not dwell on it. We have already commented on the fact that the critical density 
is now substantially larger. This is again seen explicitly in the plot of the current 
expectation value (Ji) (or (J2)), given in the left panel of figure 9. 

As far as the spatial modulation itself is concerned, we find from this analysis 
that the momentum vector kg^ for which the non-homogeneous ground state has 
minimal energy is now somewhat higher than in the Maxwell case. This is depicted 
in the right panel of figure 9. Nevertheless, it is still practically independent of the 
density. 

Numerically, the instability at pcrit ~ 3.46 corresponds, using Mj^^ = 949 MeV 
and A = 16.6 [6], to a physical density of pcrit = ~ 0.79 A^c fm~'^ when Nf = 2 
(this is close to the numerical result of [5] up to somewhat puzzling factors of two) . 

5. Discussion and open questions 

We have analysed the effect of an axial chemical potential in the Sakai-Sugimoto 
model in the low-temperature phase where chiral symmetry is broken. We have found 
that for sufficiently large chemical potential, a vectorial condensate forms, which is 
spatially modulated with a momentum vector which is practically independent of the 
potential or charge density. 

The results persist beyond the Maxwell truncation for the DBI action as well. 
The main difference between the two is that critical density and spatial modulation 
momentum take on larger values in the DBI case. The non-linearities thus in a 
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sense stabilise the homogeneous phase. It would be interesting to understand if 
the instabilities of the type we find persist if further non-linearities due to gravity 
back-reaction are included. 

It would also be interesting to understand if the new phase which we find occurs in 
other QCD-like theories, and in particular in QCD itself. The analysis of [12] does not 
find this phase, so it is worthwhile to analyse other (holographic and non-holographic) 
models and isolate more carefully which feature of the theory is responsible for the 
existence of the new phase. It would also be interesting to understand how robust 
is our finding that the spatial modulation is almost independent of the value of 
the charge density, especially since in the large- A^'c analysis for the quark chemical 
potential [1], the momentum is proportional to the chemical potential /cdgr ~ A^v- 

Putting together our results with those of Park and Ooguri [5] , we can construct 
a large part of the phase diagram of the Sakai-Sugimoto model in the presence of an 
axial chemical potential. We see that both in the confining and non-confining (and 
chirally symmetric) phases at large enough values of the axial chemical potential a 
second order phase transition to a new non-homogeneous phase appears. To complete 
the phase diagram it would be interesting to study the deconfinement in the presence 
of the non-homogeneous condensate, and see if and how it modifies the deconfinement 
temperature. For related work in other models see e.g. [14]. 

Finally, the chemical potential analysed in this paper can be embedded in the 
context discussed in [8]. The question then arises how the homogeneous (but non- 
isotropic) condensate found there competes with the non-homogeneous condensate 
found in the present paper. We will return to this question in a future publication. 
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6. Appendix: The Hamiltonian 



6.1 Maxwell truncation 



Here we present the detailed formulas for the Hamiltonian associated with the Chern- 
Simons action, taking special care about the surface terms which are present for the 
D8-D8 system. 

Our starting point is the Lagrangian density of the Maxwell-Chern-Simons action 
of (3.1), where we split the space and time indices. 



L = —K 



J' Qa + T:-' ah 



-9 



(6.1; 



The conjugate momentum associated to doAa is given by 



d{doAa) 

so that the Hamiltonian takes the form 



Oabcd A -p 
t •r\b-' cd 1 



-9 



(6.2) 



M-CS 



K / d xdz 



-9 



T-Oa T- 
-7" 7"( 



Oa 



ab 



da [^W'Ao] 



-2 



da 



-9J" 



aO\ 



3a 



Oabcd 



ab'' cd 



(6.3) 



Ao 



The first two terms in the above Hamiltonian belong to the bulk part ifbuik; the 
third term is a surface term, i.e. a boundary Hamiltonian ifbdy, while the terms in 
the last line vanish, since it is just the zeroth component of the Maxwell-Chern- 
Simons equations of motion. 

Using the Maxwell-CS equations of motion and the ansatz (3.10), the canonical 
momentum (6.2) reduces to 



H* = , H' 



2k 



9 L 



(6.4) 



and the Hamiltonian becomes Hm-cs = -f^buik + -f^bdy with 



bdy 



2KVr 



-AU p Pa 



(6.5) 
(6.6) 



where "H = M^j^YxX^Nc/lSn"^), and we have expressed the result in terms of the 
dimensionless variables defined in (3.14). We hence see that for our configuration, 
when there is non-vanishing value of chemical potential, the boundary Hamiltonian 
i^bdy is nonzero. 
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By adding to the action an extra surface term Sbdy which for our configuration 
takes the form 

Sbdy = -'^i^P j d^xdz(9^/, (6.7) 

one can remove the surface contribution to the equation of motion (3.3) which origi- 
nates from the variation of the action (3.1). This addition also simultaneously cancels 
the boundary contribution to the Hamiltonian -ffbdy When comparing free energies 
of the various configurations in this paper, we therefore just need to compute the 
bulk contribution i/buik to the Hamiltonian. 

6.2 DBI-CS Hamiltonian 

The DBl-CS Lagrangian can be written as 



C 



a 



-6 mn-F t 



The conjugate momentum for the gauge field is given by 



I'' mn-J~ pq ; 



^13 



so the corresponding Hamiltonian is 



^/Z;^£;(o,a) + 



a 



Oabcd A X" 



(6.^ 



(6.9) 



DBI-CS 



j d^xdz 17^^ 



1 + ^E^^'^^J^Oa 



+ da [^W'Ao] 



da 



7/3 



A 



(6.10) 



Note that the last term vanishes because it is just the 0th component of the DBI-CS 
equations. 

Using the DBI-CS equations and the ansatz (3.10), the conjugate momentum 
reduces to 



H* = 



H" 



1 



-9 



(— p -|- akh^) , 



(6.11) 



where p = 2Kp. This is the same as the result (6.4) found for the Maxwell truncation. 
The boundary term in the action variation takes the form 

5^bdy = P Jd^xdzdjf, (6.12) 

so in order to obtain a stationary action we need to add the boundary term 

S = -p jd'^xdzdj. (6.13) 

This term cancels the boundary term in the Hamiltonian. The bulk term takes the 
form 



with Q{z) given by (4.11). 



(6.14) 
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